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We present a study of the decays B + — » p + 7, -B — > p°7, and B° — > 1^7. The analysis is based on 
data containing 347 million BB events recorded with the BABAR detector at the PEP-II asymmetric 
B factory. We measure the branching fractions B{B + -> p+7) = (1.10±g-_|| ± 0.09) x 10~ 6 and 
B(B° -> p° 7 ) = (0.79+^^=^.06) xlO -6 , and set a 90% C.L. upper limit B(B° -* W7) < 0.78 xlO -6 . 
We also measure the isospin- averaged branching fraction B[B —> (p/uj)j] = (1.25+.Q 'jj ±0. 09) x 10~ 6 , 
from which we determine |Vtd/Vt s | = O.2OO+0020 ±0.015, where the first uncertainty is experimental 
and the second is theoretical. 

PACS numbers: 12.15.Hh, 14.40.Nd 



In the Standard Model, the decays B + — ► p+7, 
B° — > p°7, and B° — > ^7[1] arise mainly from 
b — > c?7 penguin diagrams containing a virtual top 
quark in the loop. By relating the three individ- 
ual decay rates by isospin symmetry and using the 
measured ratio between the charged and neutral 
B meson lifetimes tb+ /t^o , an isospin-averaged 
branching fraction is defined: B[B — > (p/uj)j] = 

\ {B(B+ -> p+7) + %± [B(B° - p° 7 ) + B(B° - «7)] } 

Recent calculations predict — > (,0/^)7] to be in 

the range of (0.9-1.8) x 1CT 6 [2, where most of the 
uncertainty is due to the calculation of the form factor. 
These predictions could be modified by processes beyond 
the Standard Model (3|. 

While the exclusive decay rates have a large uncer- 
tainty due to non-perturbative long-distance QCD ef- 
fects, some of this uncertainty cancels in the ratio of 
B[B — > (p/oj)^] to B — > K*-f branching fractions. Since 
the dominant diagram involves a virtual top quark, 
this ratio is related to the ratio of Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements |Vt<j/Vta| [l- HI via 



B[B -> (pA^yj 
B(B -► K*i) 



V, 



Id 



Vu 




The coefficient £ is the ratio of the form factors for the 
decays B — > ,07 and £? — > If *7 and Ai? accounts for dif- 
ferent dynamics in the decay (e.g. annihilation diagrams 
can contribute to B + — > p+7). Physics beyond the Stan- 
dard Model could affect these decays, creating inconsis- 
tencies between the measurement of |Vtd/Vt s | obtained 
from this analysis and that obtained from the ratio of 
B° and B® mixing frequencies 0]. 

Previous searches by BABAR Q and CLEO § found 
no evidence for the decays B — > pry and B — > ory. An 
observation of the decay B° — > p°7 was recently reported 
by the Belle collaboration 0]. This letter reports on a 
study of the decays B + — > p+7, _B° — > p°7, and £>° — > 
W7 based on a data sample containing 347 million BB 
events, corresponding to an integrated luminosity of 316 
fb~\ collected with the BABAR detector [Io| at the PEP- 
II asymmetric-energy e+e~ storage ring. These results 
supersede the previous BABAR measurements Q ■ 

The decays B — > p7 and B — > ujj are reconstructed 
by combining a high-energy photon with a vector me- 



son reconstructed in the decay modes p° — > ir + ir (B ~ 
100%), p+ -> tt+tt (£ ~ 100%), and u -> tt+tt-tt 
(6= [89.1 ±0.7]%) [ll|. 

The dominant source of background is continuum 
events (e+e~ — > qq, with g = u,d,s,c) that contain 
a high-energy photon from ir° or r\ decays. Other 
backgrounds include photons from initial-state radiation 
(ISR) processes, decays of B — > if*7 (if* — > Kir), decays 
of B — ► (p/ui)tt° or i? — > (p/ui)rj and combinatorial back- 
ground from higher- multiplicity 6 — > S7 decays. For each 
signal decay mode, selection requirements have been op- 
timized for maximum statistical sensitivity with assumed 
signal branching fractions of 1.0 x 10~ 6 and 0.5 x 10~ 6 
for the charged and neutral modes, respectively. 

The photon from a signal B decay is identified as 
a well-isolated energy deposit in the electromagnetic 
calorimeter with energy 1.5 < E* < 3.5 GeV in the cen- 
ter of mass (CM) frame. The energy deposit must not 
be associated with any charged track and must meet sev- 
eral other requirements designed to eliminate background 
from hadronic showers and charged particles In or- 
der to veto photons from ir° and 77 decays, we associate 
each high-energy photon candidate 7 with each of the 
other photons 7' in the event. We reject the candidates 
that are consistent with originating from ir° or rj decays 
based on a likelihood ratio constructed from the energy 
of the second photon 7' and the invariant mass of the 
pair m 77 ' . We also combine the high-energy photon can- 
didate with photon conversions to e+e~ pairs, and reject 
the photon if the invariant mass is consistent with a 7r° 
or 77. 

Charged-pion candidates are selected from well- 
reconstructed tracks with a minimum momentum trans- 
verse to the beam direction of 100 MeV/c. A stringent 
selection algorithm is applied to reduce background 
from charged kaons produced in b — > sj decays. The al- 
gorithm combines the information provided by the ring- 
imaging Cherenkov detector with the measurement of en- 
ergy loss in the tracking system. 

Photon candidates with energy greater than 50MeV 
in the laboratory frame are combined into pairs to form 
7T° candidates. For B° — > u>7 (B + — > p+7) decays, the 
invariant mass of the pair is required to satisfy 122 < 
m 77 < 150 MeV/c 2 (117 < m 77 < 148 MeV/c 2 ). We 
also require that the cosine of the opening angle between 
the daughter photons in the laboratory frame be greater 



5 



than 0.413 (0.789). 

The identified pions are combined into vector meson 
candidates by requiring 633 < m^+„- < 957 MeV/c 2 , 
636 < m I+I o < 932 MeV/c 2 , and 764 < m^+^-^o < 
795 MeV/c 2 for p°, p + , and oj, respectively. The charged 
pion pairs must originate from a common vertex. The 
separation along the beam axis between this vertex and 
the one obtained by combining the other charged parti- 
cles in the event is required to be less than 4 mm and to 
be measured with a precision better than 0.4 mm. 

The photon and p/uj candidates are combined to form 
the B meson candidates. We define AE = E* B — £^ eam , 
where E* B is the CM energy of the B meson candidate and 
^"beam * s the CM beam energy. We also define the beam- 
energy-substituted mass toes = \/ -S'bcam ~ Pb^ ' wn ere 
p B is the CM momentum of the B candidate. Signal 
events are expected to have a AE distribution centered 
near zero with a resolution of about 50 MeV, and an toes 
distribution centered at the mass of the B meson, with 
a resolution of 3 MeV/c 2 . We consider candidates in the 
ranges -0.3 < AE < 0.3 GeV and m ES > 5.22 GeV/c 2 , 
which include sidebands that allow the combinatorial 
background yields to be extracted from a fit to the data. 

In signal events the vector meson is transversely polar- 
ized, while in background events from B — * p(ir° jr\) and 
B — > ti>(7r°/r/) it is longitudinally polarized. To reject this 
background, we calculate the vector meson helicity angle, 
9 H , and require | cos9 H \ < 0.75. The helicity angle is de- 
fined as the angle between the B momentum vector and 
the 7T~ track calculated in the p rest frame in the case of 
a p meson, or the angle between the B momentum vector 
and the normal to the u> decay plane for an uj meson. 

Contributions from continuum background processes 
are reduced by considering only events for which the ra- 
il^ of second-to-zeroth Fox- Wolfram moments 13 1 
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calculated using the momenta of all charged and neu- 
tral particles in the event is less than 0.7. A neural 
network combining the variables described below further 
suppresses continuum background. The quantity R' 2 , de- 
fined as i?2 in the frame recoiling against the photon 
momentum, is used to reject ISR events. To discrim- 
inate between the jet-like continuum background and 
the more spherically-symmetric signal events, we com- 
pute the angle between the photon and the thrust axis 
of the rest of the event (ROE) in the CM frame. The 
ROE is defined by all charged tracks and neutral energy 
deposits in the calorimeter that are not used to recon- 
struct the B candidate. We also calculate the moments 
Li = J2j Pj ■ I cos 0j I 7 J2j P*j i where p* and 9* are the mo- 
mentum and angle with respect to an axis, respectively, 
for each particle j in the ROE. We use Li, Li, and L3 
with respect to the thrust axis of the ROE, as well as 
with respect to the photon direction. In addition, we 
calculate the B meson production angle 9* B with respect 
to the beam axis in the CM frame. Differences in lepton 
and kaon production between background and B decays 



are exploited by using flavor-tagging variables [14]. The 
significance of the separation along the beam axis of the 
B meson candidate and ROE vertices is included as well. 
The purity of the selected sample is enhanced by a cut on 
the output of the neural network that retains 63%, 74%, 
and 71% of the signal events in the modes B + — > p + -f, 
B° — > p°j, and B° — ► cwy, respectively. 

The expected average candidate multiplicity in the se- 
lected signal events is 1.01 for B° — > p 0r y and 1.07 for 
both B + — ► p + 7 and B° — > wy. In events with multiple 
candidates, the one with the reconstructed vector meson 
mass closest to the nominal mass is retained. This crite- 
ria was chosen because the mass of the vector meson was 
found to be uncorrelated with the variables used in the 
fit. Applying all the selection criteria described above, 
we find efficiencies [H] of 11.0% for B+ -> p+j, 14.1% 
for B° -> p°7, and 7.9% for B° -> luj. 

The signal content of the data is determined by a mul- 
tidimensional unbinned maximum likelihood fit, which is 
constructed individually for each of the three signal decay 
modes. All fits use AE, toes, cos9h, and the neural net- 
work output N. In order to facilitate the parametrization 
of the probability density function (PDF) used in the fit, 
the transformation AW = tanh -1 (ci • N — C2), in which 
the Ci are are mode- dependent constants, is made. For 
decays B° — > loj (lo — > 7r + 7r~7r°), the cosine of the angle 
between the 7r + and 7r° momenta in the 7r + 7r~ rest frame 
(Dalitz angle) is added as a fifth observable. 

In the fit we consider several hypotheses for the origin 
of the events: signal, continuum background, B — > K*j 
decays, and other B backgrounds. The likelihood func- 
tion for a signal mode k (= p +r y, p°7, 0*7) is defined as 



A'h 



exp 



E 




, (2) 



where A^h yp is the number of event hypotheses, n, is the 
yield of each hypothesis, and Nk is the number of can- 
didate events observed in data. Since the correlations 
among the obscrvablcs are found to be small in simulated 
event samples, we define the PDF Ti{xj;ai) for the i-th 
event hypothesis as the product of individual PDFs for 
each fit observable Xj given the set of parameters cti. 

Each PDF is determined from a one-dimensional fit to 
a dedicated sample of simulated events. The AE PDF 
is corrected for the observed difference between data and 
simulation by using samples of B — > if* 7 decays. All 
continuum background parameters float freely in the fits, 
while the shapes of the signal and B background distribu- 
tions are fixed according to the Monte Carlo simulation. 
The signal toes spectra are described by Crystal Ball 
functions the angular distributions are modeled by 
polynomials, and the distributions of AE and AW are 



parametrized as f(x) — exp 



where 



p is the peak position of the distribution, <Jl,r are the 
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FIG. 1: AE and rn.ES projections of the fits for the decay 
modes B + -> p + j (top), 5° -» p°j (middle), and B° -> lo^ 
(bottom). In each plot, the signal fraction is enhanced by 
selections on the other fit variables. The points are data, the 
solid line is the total of all contributions, and the long-dashed 
(dashed-dotted) line is background-only (signal-only). 



TABLE I: The signal yield (n s i g ), significance (£) in stan- 
dard deviations including systematic errors, efficiency (e), and 
branching fraction (B) for each mode. The errors on n s j g are 
statistical only, while for the branching fraction the first error 
is statistical and the second systematic. 



Mode n s i g 


E 


e(%) 


B(10" 


-6 




B+ p+ 7 42.0t£? 


3.8a 


11.0 


1.10 


+0.37 
-0.33 


± 


0.09 


B° - A 38.7^° 8 6 


4.9a 


14.1 


0.79 


+0.22 
-0.20 


± 


0.06 




2.2a 


7.9 


0.40 


+0.24 
-0.20 


± 


0.05 


B -> {p/u;)j 


6.4a 




1.25 


+0.25 

-0.24 


± 


0.09 


B -* p-y 


6.0a 




1.36 


+0.29 
-0.27 


± 


0.10 


TABLE II: Fractional systematic errors 


(in %) of 


the mea- 


sured branching fractions 














Source of error 


p+7 


A 


OJ'f 


PI 




(p/w)7 


Tracking efficiency 


f.O 


2.0 


2.0 


1.4 




1.5 


Particle identification 


2.0 


4.0 


2.0 


2.9 




2.7 


Photon selection 


1.9 


2.6 


1.7 


2.2 




2.1 


7T° reconstruction 


3.0 




3.0 


1.9 




2.5 


7T° and r\ veto 


2.8 


2.8 


2.8 


2.8 




2.8 


MM efficiency 


1.0 


1.0 


1.0 


1.0 




1.0 


MM shape 


0.4 


0.3 


2.3 


0.4 




0.7 


Signal PDF shapes 


4.8 


3.3 


2.4 


3.1 




2.6 


B background PDFs 


3.9 


2.9 


9.7 


3.2 




3.1 


BB sample size 


1.1 


1.1 


1.1 


1.1 




1.1 


B(LO -> TT + 1V-TT ) 






0.8 






0.1 


Sum in quadrature 


8.1 


7.4 


11.6 


7.0 




6.9 



widths on the left and right of the peak, and (Xl,r are a 
measure of the tails on the left and right of the peak, re- 
spectively. Various functional forms are used to describe 
the continuum and B background components. 

We measure the signal yield n S i g by maximizing the 
likelihood function in Eq. ([2]). In the fit, the contin- 
uum background yield is allowed to float, as is the over- 
all yield of the B background, with the exception of the 
B+ -► K*+~f {K*+ -> K+tt°) yield in the B+ p+j 
mode, which is fixed. The relative yields among the dif- 
ferent B backgrounds are fixed to the values obtained us- 
ing known branching fractions ll| and selection efficien- 
cies determined from simulated events. Figure [1] shows 
the data points and the projections of the fit results for 
AE and toes separately for each decay mode. The signal 
yields are reported in Table |TJ The significance is com- 
puted as \/2A log£, where Alog£ is the log-likelihood 
difference between the best fit and the null-signal hypoth- 
esis. 

Table (Til gives an overview of the contributions to the 
systematic uncertainties. These are associated with the 
signal reconstruction efficiency and the modeling of sig- 
nal and BB background in the Monte Carlo simulation. 



The latter contributes to the uncertainties on the signal 
yields. The systematic error affecting the signal effi- 
ciency includes uncertainties on tracking, particle iden- 
tification, 7 and 7r° reconstruction, ir°/r] veto and the 
neural network selection. The uncertainties on the ir° /i] 
veto and neural network selection are determined from 
a control sample of B — > Dir decays, with D — > Kit or 
D — > Kirn. To estimate the uncertainty related to the 
modeling of the signal and B background in the Monte 
Carlo, we vary the parameters of the PDFs that are fixed 
in the fit within their errors. The uncertainty related 
to the choice of a specific functional form for the shape 
of the MM distribution is evaluated by using a binned 
PDF as an alternative description. All relative and ab- 
solute normalizations of B background components that 
are fixed in the fit are varied within their errors. For all 
these variations, the corresponding change in the fitted 
signal yield is taken as a systematic uncertainty. 

The branching fractions are calculated from the mea- 
sured signal yields assuming B(T(4S) — ► B°B°) = 
B{T(AS) -> B+B~) = 0.5. The results are listed in Ta- 
blelU For B° — > uij, we also compute the 90% confidence 
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level (C.L.) upper limit B(B° -> ujj) < 0.78 x 10~ 6 using 
a Bayesian technique. We determine the branching frac- 
tion upper limit Bi such that J®' CdB/ J °° LdB = 0.90, 
assuming a flat prior in the branching fraction and taking 
into account the systematic uncertainty. 

We test the hypothesis of isospin symmetry by measur- 
ing the quantity T(B+ -> p+-f)/[2T{B° -> p° 7 )] - 1 = 
— 0.35±0.27. The result is consistent with the theoretical 
expectation 0. 

The isospin-averaged branching fraction is extracted 
from a simultaneous fit to the three decay modes: 

B[B -> (p/w)7] = (1.25±g;|| ± 0.09) x 10~ 6 . (3) 

In the fit we impose the isospin constraints on the widths 
of the decay modes: r B ^ p + 7 = 2T B ^ p a 1 = ZTb^u-i- 
Our measurements of the individual branching fractions 
are consistent with this hypothesis with a x 2 of 1.8 for 
2 degrees of freedom. The significance of the signal is 
6.4cr, including systematic uncertainties. This result is 
consistent with the measurement from Belle [9j. If we 
exclude the B° —>■ ujj mode from the simultaneous fit, 
we obtain B(B -> fry) = (1.36^2? ± 0.10) x 10" 6 . Using 
the world average value of B(B — ► K*j) [11| . we calculate 
B[B -► {p/uj)i\/B{B -> AT* 7 ) = 0.030 ± 0.006. This 
result can be used to calculate the ratio |Vtd/V ts | 0, 
Hjj ]. Following [lj}, we obtain 

\Vtd/V u \ = 0.200l°;°2J ± 0.015, (4) 

where the first error is experimental and the second is 
theoretical. This result is consistent with the measure- 
ment of this ratio from the study of B° and B® mixing Q . 

In conclusion, we have measured the branching 
fractions of B{B+ -> p+ 7 ) = (1.10±g-Jl ± 0.09) x 10~ 6 
and B(B° -> p° 7 ) = (0.79±gjg ± 0.06) x 10~ 6 , and 
set a 90% C.L. upper limit on the B° — ► w 7 branch- 
ing fraction of 6(5° -> w 7 ) < 0.78 x 10~ 6 . The 
isospin-averaged branching fraction B[B — > (p/u;) 7 ] = 
(1.25^Q 24±0.09) x 10 -6 is the most precise measurement 
of this quantity to date. This measurement is used to 
extract \V td /V ts \ = 0.200±°;° 2 J ± 0.015. 
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